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We study a distinct supersymmetric signal of multi-photons in association with jets and missing 
transverse energy. At least one of these photons has the origin in displaced vertex, thus delayed 
and non-pointing. We consider a supersymmetric scenario in which the gravitino is the lightest 
supersymmetric particle (LSP) (with a mass ~ 1 keV) and the lightest neutralino is the next- 
to-lightest supersymmetric particle (NLSP). The NLSP decays dominantly into a photon and a 
gravitino within the detector with a decay length ranging from ct^ ~ 50-100 cm. In addition, we 
assume that the second lightest neutralino and the lightest neutralino are nearly degenerate and 
this leads to a prompt radiative decay of the next-to-lightest neutralino into a photon and a lightest 
neutralino with a large branching ratio. Such degenerate neutralinos can be realised in various 
representations of the SU(5), SO(10), and E(6) Grand Unified Theories (GUTs). The non-pointing 
photons can be reconstructed at the electromagnetic calorimeter of the ATLAS inner-detector, which 
have been designed with good timing and directional resolution. We find that with a centre-of-mass 
energy E cm = 14 TeV at an integrated luminosity of 100 fb~ l one may see evidence of hundreds 
of tri-photon events and a few four-photons events at the LHC, in addition to several thousands 
di-photon events. We also predict the event rates even at the early phase of LHC run. 

PACS numbers: 12.60.Jv, 14.80.Ly, 14.80. Nb, 11.30.Pb 



I. INTRODUCTION 

In this era of the Large Hadron Collider (LHC) the 
TeV scale physics is expected to be probed. Supersym- 
metric Standard Model (SSM) is one of the most inter- 
esting and attractive candidate for physics beyond the 
Standard Model (SM). It offers a possibility of gauge 
coupling unification and dark matter candidate, and also 
solves the gauge hierarchy problem. Once supersymme- 
try (SUSY) is realized as a local symmetry [1], it predicts 
the existence of the gravitino G as the spin-3/2 super- 
partner of the graviton. Supersymmetry breaking leads 
to a non-zero mass of the gravitino through the super- 
Higgs mechanism, in which the gravitino "eats up" the 
spin- 1/2 goldstino associated with spontaneously broken 
local supersymmetry [2-5] . The mass m,Q of the gravitino 
is governed by the scale of SUSY breaking and can range 
from as low as eV scale to as high as 100 TeV scale [6-13]. 
In this work we choose a phcnomenological supersymmet- 
ric scenario in which gravitino is the lightest supersym- 
metric particle (LSP) with a mass ijiq ~ 1 keV and look 
at the collider signatures of such a scenario at the LHC. 
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Such light gravitinos also have implications in cosmology. 
First of all, one should note that the dark matter relic 
density is presently known to be VtoMh 2 ~ 0.11 [14]. 
In addition, constraints on structure formation require 
that the bulk of the dark matter should be cold or warm 
[15]. For a gravitino with a mass ~ 1 keV, nonstan- 
dard cosmology and a nonstandard gravitino production 
mechanism are required to satisfy small-scale-structure 
constraints and to avoid overclosure [16]. One might also 
need some other dark matter particle. An example of a 
nonstandard early-Universe physics is to consider a low- 
reheating temperature [17, 18]. In Ref.[18] a low-reheat 
scenario has been proposed in which a gravitino of mass 
m,Q = 1-15 keV can have the right abundance to be the 
warm dark matter. 

The interactions of the gravitino are suppressed by the 
reduced Planck Scale M P = 2.4 x 10 18 GeV and a light 
gravitino interacts more strongly than a heavy gravitino. 
Light gravitinos are primarily produced at colliders in 
the decays of the NLSP. In our scenario the lightest neu- 
tralino (xi)- which is predominantly a bino, is the NLSP 
and it decays dominantly into a photon and a gravitino. 
These photons are delayed and non-pointing as they are 
not pointing to the interaction vertex where the NLSP 
is produced. Along with this we also look into the ra- 
diative decay of the second-lightest neutralino (x°) i.e., 
X2 X?7> where the emitted photons are prompt. Thus 
our main goal in this paper is to study the spectacular 
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multi-photon events at the LHC where there is a mixture 
of prompt photons and non-pointing photons in the final 
states. In order to have a large branching ratio of the 
decay x.2 —> X?7j we choose a framework where the U(l) 
and SU(2) gaugino soft SUSY breaking mass parame- 
ters Mi and M2 , respectively are very close and result in 
nearly mass degenerate x.2 anc ^ Xv 

In a minimal supergravity like framework (mSUGRA) 
the gaugino masses are unified at the high scale (unifica- 
tion scale) . When they run down to electroweak symme- 
try breaking scale (EWSB) the gaugino mass ratio gets 
modified through renormalisation group effects (RGEs). 
At the EWSB scale the approximate ratio of the gaugino 
masses are given as M\ : M 2 : M 3 ~ 1 : 2 : 6, where 
M3 is the SU(3) gaugino soft SUSY breaking mass pa- 
rameter and M\ , M2 have been defined in the previous 
paragraph. So it is very clear from the above ratios that 
in a mSUGRA scenario it is almost impossible to have 
nearly degenerate neutralions at the EWSB scale. But if 
the gauginos masses are non-universal at the high scale 
with M\ > M2 then the RGEs compensate for M2 and 
one can have nearly degenerate gauginos at the EWSB 
scale. In this paper we point out a few grand unified 
gauge symmetry breaking patterns where this feature can 
be grabbed. 

Light gravitino and its collider signatures have been 
studied extensively in various context [16, 19-44] and 
mostly in connection with gauge mediated supersym- 
metry breaking (GMSB) [8-10]. Signatures involving 
photons are characteristics of scenarios with neutralino- 
NLSP. In most of the cases studied so far the lightest 
neutralino is predominantly a bino and the second light- 
est neutralino is dominated by its wino component with a 
large mass splitting between them. However, as empha- 
sized earlier, we will consider a scenario where the light- 
est and the second lightest neutralino are approximately 
degenerate in mass. This will lead to multi-photon signa- 
tures at the LHC for a 1 keV gravitino, where in the final 
states we can have combinations of prompt and delayed 
photons. This is a spectacular signal free from Standard 
Model backgrounds and has not been studied earlier. The 
signature of two non-pointing photons is very much dis- 
tinct and clean with a large event rate at the LHC. We 
discuss the di-, tri-, and four-photon signals at 14 TeV 
center-of-mass (CM) energy with 100 fb _1 integrated lu- 
minosity. We find it very hard to get any significant 
event rates for four photons at 7 TeV CM energy with 
3 fb _1 , and it is not a surprise. Let us note in pass- 
ing that triphoton signatures of Randall-Sundrum model 
have been studied recently in Ref. [45] . 

We discuss the pr distributions of multi photons for 
different suggested benchmark points (BP). We use the 
decay kinematics of the neutralino with a sufficiently long 
lifetime. Schematic diagram of a neutralino decaying into 
a gravitino and a photon in the ATLAS detector is shown 
[35] in Fig. 1. If the decay length of the Xi IS comparable 
to the size of the ATLAS inner-detector [35, 46], high- 
Pt photons could enter the calorimeter at angles (?7 7 ) 



deviating significantly from the nominal angle from the 
interaction point to the calorimeter cell (rji ) . 
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FIG. 1: Decay kinematics of the NLSP (the lightest neu- 
tralino) in the ATLAS detector [35, 46]. 

The plan of the paper is as follows. We discuss the 
gravitino production from Xi decay in Sec. II. In Sec. 
Ill we discuss how nearly degenerate gaugino masses en- 
hance the branching ratio of the radiative decay of the 
next-to-lightcst neutralino and suggest the possible high 
scale scenarios from where this degeneracy condition can 
be achieved. We suggest four benchmark points for our 
numerical analysis and these are discussed in Sec. IV. 
In Sec. V we start with our goal for collider simula- 
tion and discuss the multi-photons+jeis associated with 
missing transverse energy (MET) as a potential signal at 
the LHC. Sec. VI contains the results of our numerical 
analysis and the conclusion is provided in Sec. VII. 



II. GRAVITINO PRODUCTION FROM 
NEUTRALINO DECAY 

As discussed earlier, the gravitino gets a mass by the 
super-Higgs mechanism. The mass of the gravitino is re- 
lated to the fundamental supersymmetry-breaking scale 
y/F, as 



F 



V3M F 



240 eV 



F 



10 3 TeV 



(1) 



The weak-scale gravitino has a very feeble interaction and 
thus it is usually hard to find its signatures in collider ex- 
periments. However, once SUSY is broken spontaneously, 
the extremely weak gravitino interactions are enhanced 
at energy scales much larger than the gravitino mass tuq. 
This is because in the high energy limit the gravitino has 
the same interaction as the goldstino and the couplings 
of the goldstino are proportional to 1 /F [47-49] . Hence 
the decays of heavier sparticles to gravitinos are faster 
for light gravitinos. The partial decay widths of Xi to G 
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are given as [22, 32, 50]: 
fci 



r(x? 



r(x? -> </>g) 

where, 
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(3) 
(4) 



| Ni i cos 6 W + N 12 sin 6, w \ 2 , 
liVusin^ + iViacos^l 2 , 
|iVi3 cos Bp — Nu sin#g| 2 , 
|JVi3 sin a — N14 cos a\ 2 , 
\Ni3 cos a + N14 sina| 2 , 
IJV13 sin,3 + iVi 4 cos/?| 2 . 



(5) 



Here Nij are the neutralino mixing matrices, 9 W is the 
weak mixing angle, a is the Higgs (<I> = h°, H°, A ) mix- 
ing angle and tan j3 is the ratio of the vacuum expectation 
values of the two Higgs doublets Hi and H 2 in the SSM. 
From the above expressions it is clear that for a bino-likc 
NLSP Nncos6 w is much larger than .ZV12 sin^. The 
decay modes into the photon dominates over Z and (f> 
channels as the later two have phase-space suppressions. 
Assuming that the decay widths in Z and 4> channels are 
negligible, the decay length of the lightest neutralino is 
given by 




100 GcV 



F 



100 TeV 



x 10~ 2 cm. (6) 



For a pure bino-likc lightest neutralino and nig = 1 kcV, 
we get a decay length ct^ « 70 cm. 



III. RADIATIVE DECAY OF NEUTRALINO: 



X°2 



xh 



The radiative decay of second lightest neutralino em- 
anates at the one-loop level and decay width is given as 
[511: 
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where g^o^o^ oc eg 2 /16ir 2 is an effective coupling. This 
radiative decay is enhanced [51, 52] by a kinematic factor 
when x1 and X2 are nearly degenerate as in this regime 
of parameter space three body decays are suppressed by 
a factor e 5 , where e = (1 — m^o/m^o). It is being noted 
in [52-54] that the decay branching ratio of \2 ~~ ^ xh 
is much larger for large tan j3 and /1 > Mi , M 2 with 



~ Mitan/3/2. In the general MSSM scenario the 
radiative decay branching ratio can reach nearly 100% 
[54] for [Mi|, M 2 < 1000 GeV, with |Mi|, M 2 < 
and I Mi | ~ M 2 . We have calculated the branching ratio 
of the radiative decay of the second lightest neutralino 

using SDECAY version 1.3b [55]. 



A. Radiative decay with non-universal gaugino 
masses 

In [54] , the enhancement conditions in the radiative de- 
cay branching fractions arc justified for the minimal su- 
pergravity(mSUGRA) models with non-universal gaug- 
ino masses. The part of the N=l supergravity La- 
grangian (the part that contains only the real part of 
the left-chiral superficlds $j) containing the kinetic en- 
ergy and the mass terms for the gauginos and the gauge 
bosons can be written as: 



l C = 



- \Ref a ^){-l/2\ a If)\P) - \Ref a ^)F^F^ 

+ ie- G / 2 G l ((G- 1 )i)[a/^(0*)/9^']A Q A /3 + h.c 

(8) 

where G % — dG/d<j>i and (G" 1 )* is the inverse matrix 
of G^i = dG/d(f>* l dcf)j , A Q is the gaugino field, and cf> is 
the scalar component of the chiral supcrficld $ and F^ v 
is the unified gauge kinetic term. The F-component of 
the symmetry breaking scalar field $ generates gaugino 
masses with a consistent SUSY breaking with non-zero 
vacuum expectation value (vev) of the chosen F, where 



F j = ie- G / 2 [G*((G-^')]. 



(9) 



The $ J 's can be a set of GUT singlet supermultiplets 
<& s , which are part of the hidden sector, or a set of non- 
singlet ones $ N , fields associated with the spontaneous 
breakdown of the GUT group to SI/ (3) ® SU(2) <g> U{1). 
The non-trivial gauge kinetic function f a f}{&) can be 
expanded in terms of the non-singlet components of the 
chiral superficlds in the following way 



N 



+ 0(—) 2 

Mp + { Mp> ' 



a/3 



(10) 



where /o and £ N are functions of chiral singlet super- 
fields, essentially determining the strength of the inter- 
action and Mp is the reduced Planck mass= Mpi / \/8n. 
The contribution to the gauge kinetic function from <f> N 
has to come through symmetric products of the adjoint 
representation of the associated GUT group, since f a p 
has such transformation property for the sake of gauge 
invariance. The non-universal gaugino masses are calcu- 
lated for SU(5), 5*0(10) and E(6) grand unified gauge 
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Representations 


Mi : M2 : M 3 
fat AfniT^ 


Mi : M 2 : Af 3 
fat Mt\ 


75 C 577(5) 
210, 770 C SO(10) 


-5:3:1 


-5:6:6 


2430 C E(6) 


-1=1:1 

5 . 3 . 1 

2 ■ 2 • 


-1.8 : 2 : 6 
2.5 : -3 : 6 



TABLE I: Ratios of gaugino masses for F '-terms in represen- 
tations of 5(7(5), 50(10) and E(6) leading to nearly degen- 
erate gauginos at low scale. 



groups in [56]. The results for the ratios of gaugino 
masses are given in Table I. We have tabulated here only 
the cases where Mi and M2 are nearly degenerate at the 
EWSB scale with Mi < M<z. This fits in our scenario. 



IV. BENCHMARK POINTS 

In this section we present four benchmark points (see 
Table II) we have worked with to demonstrate that the 
nearly degenerate Mi and M2 at the EWSB scale can 
lead to radiative decay of the second lightest neutralino. 
In addition to this we have also shown that if one has 
Mi < M2 and the gravitino in the bottom of the spec- 
trum, then the lightest neutralino has a sizeable branch- 
ing fraction to decay into a photon and gravitino. This 
leads to multi-photon signatures in collider experiment. 
The spectrum has been generated using the SuSpect ver- 
sion 2.41 [57] with all the input parameters specified at 
the electroweak scale. The gravitino mass is taken to be 
~ 1 keV which is necessary for the fact that the light- 
est neutralino decays within the detector with a decay 
length ct^o ~ 50-100 cm. The radiative decay of the 
X2 an d decay of X? have been calculated using SDECAY 
version 1.3b [55]. The benchmark points we have worked 
with are consistent with all the low energy constraints 
like muon (g — 2)^, b — > 57 and the LEP limit on the 
lightest Higgs boson mass and other charged particles 
masses [58, 59]. 

Throughout all of our benchmark points we have kept 
the value of m^o and m^o nearly the same with different 

A2 Xl J 

choices of fx, tan p, squarks, gluino and slepton masses. 
The high value of fi is important for enhancement of the 
radiative decay branching fraction of X2 m to a x?7 pair. 
We have worked with nig starting from as low as 413 GeV 
to 740 GeV. We set A t = A T = A b = A = -1000 GeV. 
The large value of \Aq\ is required for obtaining a large 
radiative decay branching ratio (BR) of x.2- For Aq~ 
0, the three-body-decay modes of X2 are dominant and 
the radiative decay is very much suppressed in our case. 
We have also noted that the radiative decay branching 
fraction depends less significantly on the sign of Aq. It is 
a little less for the positive value of Aq than the negative 
one keeping |^4o| same. In Table III we tabulate the decay 
branching fraction of the X2 an d Xi f° r our choice of 



input parameters. From this table one can see the effect 
of [i, tan/3, squark and slepton masses on the radiative 
decay of X2- However, the decay branching fraction of 
the lightest neutralino is determined once the mass of 
the gravitino, the lightest neutralino and the neutralino 
mixing parameters are fixed and does not depend at all 
on the choices of squarks, gluino and sleptons masses. 
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TABLE II: Proposed benchmark points (BP) for the study of 
radiative decay of X2 ond the NLSP \i- We have set Ao = 
— 1000 GeV for the third generation squarks and sleptons and 
it is zero for the rest. Masses of the particles and /j, are given 
in GeV. 
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0.30 


0.11 


0.26 


0.10 


Xi -> lG 


0.89 


0.89 


0.87 


0.88 



TABLE III: Branching fractions for the decays x% ~ * Xil an ^ 
X? — > 7G for different benchmark points. 



V. COLLIDER SIMULATION 

The Xi an d X2 are produced in cascade decays of 
squarks and gluinos accompanied by hard jets. In an 
i?-parity conserving scenario the gravitino is produced 
at the end of each cascade, which goes undetected at 
the collider detector, leading to large amount of miss- 
ing transverse energy (fir) (see, Fig. 2). Thus one can 
have multi-photon signals in association with hard jets 
and f$T- The collider simulation has been done with a 
centre of mass energy E cm =14 TcV, at an integrated lu- 
minosity of 100 fb _1 using the event generator PYTHIA 
6.4.16 [60]. A simulation for the early LHC run at E cm — 7 
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FIG. 2: Missing Energy distribution for different benchmark 
points with E crn — 1A TeV. 



TeV and integrated luminosity of 3 fb _1 has also been 
performed. We have used the parton distribution func- 
tion CTEQ5L [61] with the factorisation (hf) an d renor- 
malisation (/j,r) scale set at hr = [ip =average mass of 
the final state particles produced in the initial hard scat- 
tering. The effects of Initial and Final State Radiation 
(ISR/FSR) have also been taken into account. Below we 
mention the numerical values of various parameters used 
in our calculation [58] 

M z = 91.187 GeV, M w = 80.398 GeV, M t = 172.3 
GeV, a-^{M z ) = 127.9, a s {M z ) = 0.118, where M z , 
Mw and M t are the masses of the Z-boson, IF-boson 
and top quark, respectively. a em (Mz) and a s (Mz) are 
the electromagnetic coupling constant and strong cou- 
pling constant respectively at the scale of Mz- 
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FIG. 3: pt distributions of the prompt photon (top) and non- 
pointing photon (bottom) with E cm =lA TeV for all bench- 
mark points. 



A. Event selection criteria 

We have considered the following final states to demon- 
strate the event rates in multi-photon channels: 

• 2-/ + fi T + jets 

• 37 + fir + jets 

• 47 + pT + jets 

where at least one of these photons has the origin in 
displaced vertex due to the fact that the decay length 
of the lightest neutralino is O (50-100 cm). The photon 
out of a x° decay is soft (see Fig. 3-top) while the pr of 
the photon coming from a \i are normally hard (see Fig. 
3-bottom) as the mass difference between X2 an d Xi i s 
O(30 GeV). 

The following requirements have been implemented to 
select isolated photons: 

• We have identified photons with px more than 30 
GeV and |rj| <2.5. 



• A minimum AR separation between two photons 
has been demanded in t erms of AR > 0.2, where 
AR= v /(Ayy) 2 + (Aci) 2 . 

• A lepton-photon and jet-photon isolation of 
Ai?; 7 > 0.4 and ARjy > 0.6, respectively have 
been imposed. 

• The sum of hadronic Et deposit in a cone of AR = 
0.2 around the photon is required to be £|£t| < 
10 GeV. 

• To reduce the di-photon background from 7r° — > 27 
we have also required a photon-photon invariant 
mass cut - 20 GeV < M 77 < + 20 GeV. 

The photons have been ordered according to their 
hardness (see Figs. 4, 5 and 6) and a minimum px cut has 
been imposed on each of them depending on the various 
final states: 

• di-photon: p T ^ t > 50 GeV, pr 72 > 40 GeV 

• tri-photon: p Ti > 50 GeV, p Tl > 40 GeV , 
pr > 30 GeV 
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FIG. 4: pt distributions of di-photons for (from top to bot- 
tom) BP-1, BP-2, BP-3, and BP-4 with E cm =lA TeV. 



FIG. 5: pt distributions of tri-photons for (from top to bot- 
tom) BP-1, BP-2, BP-3, and BP-4 with E cm =U TeV. 



7 




• four-photon: p? yi 
p T -, 3 > 30 GeV, v T _ ti 



> 50 GcV, pt,. 

> 30 GeV. 



> 40 GcV 



Wc have also incorporated the probability of jet-faking 
as photon, which is taken to be 0.1% [46, 62] and an 
identification efficiency of 60% has been used for the non- 
pointing photons following [46, 63]. We have not taken 
into account the rapidity dependence of the identification 
efficiency and used a uniform efficiency for a conservative 
approach. 



VI. RESULTS 
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FIG. 6: pt distributions of four photons for (from top to 
bottom) BP-1, BP-2, BP-3, and BP-4 with £' cm =14 TeV. 



In this section we present the numerical results of our 
simulation. In Table IV we have presented the number 
of events in the multi-photon channels after applying the 
basic cuts listed in the previous section. The different 
benchmark points we have selected correspond to similar 
m^o and m^o with different values of /i, tan /3, slepton, 
and squark masses as given in Table II. The radiative 
decay branching fraction of x° ~~ ► Xi7 ( see Table III) de- 
pends on the choice of squarks and slepton masses as well 
as on the values of tan /3 and /i, which in turn affect the 
event rates in various multi-photon channels. In BP-1 
the gluino is lighter than the squarks. In this case, \2 
is produced via radiative and three-body decay of gluino 
which together has a branching fraction of more than 
50%. The left-handed squarks decay into a x^g-pair ei- 
ther directly (with a branching fraction ~ 8%) or via 
gluino decay. The right-handed squarks mainly decay 
into a gluino and a quark pair and the gluino further can 
decay into a X29 or X2<79~P arr - The situation is similar 
in BP-2 with only difference is that it has smaller ra- 
diative decay branching fraction (11%) of decaying into 
a X2 — Xi7"P arr due to small tan/3 = 10. In BP-3 the 
squarks are lighter than the gluino. In this case the gluino 
directly decays into a <j<j-pair. Therefore the production 
cross-section of x° in SUSY cascade decreases as the dom- 
inant contribution in this case comes only from the decay 
of left-handed squarks with a branching fraction ranging 
from 30%-35%. The radiative decay branching fraction 
X2 — * X?7 is slightly greater than BP-2, due to the fact 
that the squarks and slepton masses are smaller than that 
in BP-2 (see Table II), which contribute in the loop. 

Above all, due to different squarks and gluino masses at 
different benchmark points the overall SUSY production 
cross-section changes from one benchmark point to the 
other. This combined with the different decay branching 
fractions of qi, — > X2Q f° r various benchmark points, af- 
fects the production cross-section of the second lightest 
neutralino in cascade decay of squarks and gluino and 
shows up in the final event rates. 

From Table IV one can find that in the di-photon chan- 
nel one has substantial rate at E cm = 14 TeV and inte- 
grated luminosity of 100 fb" 1 . The dominant contribu- 
tion to this channel comes from the two non-pointing 
photons out of a x? Gj decay from the two cascades, 
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SIGNAL 


BP-1 


BP-2 


BP-3 


BP-4 


27 + fir + jets 


7942 


6604 


8150 


9549 


37 + p T + jets 


597 


220 


165 


162 


47 + p T + jets 


8 


2 


3 


3 



TABLE IV: Number of signal events, after applying the basic 
cuts at an integrated luminosity of 100 fb^ 1 and the center of 
mass energy of 14 Te V for all our benchmark points. 

and constitutes of more than 92% of the total di-photon 
cross-section. The decay branching fractions of Xi ~> O7 
are more than 87% for all of our benchmark points. The 
next sub-dominant contribution to it comes from one 
non-pointing photon from Xi decay and the other prompt 
photon from radiative decay of x.2- This constitutes ~ 7% 
of the total di-photon cross-section. The rest comprises 
of two prompt photons when we have radiative decay of 
X2 from both the cascade or a combination of prompt 
or non-pointing photon with the ISR/FSR photon, the 
fraction of which is rather small. We have also presented 
the di-photon rates even at the early phase of LHC run 
with E cm = 7 TeV at an integrated luminosity of 3 fb _1 
(see Table V). In BP-1, the di-photon rate is larger than 
BP-2 with nearly identical spectrum. This attributes to 
the fact that the radiative decay branching fraction of X2 
at BP-2 is one-third of that in BP-1. 



SIGNAL 


BP-1 


BP-2 


BP-3 


BP-4 


27 + p T + jets 


20 


17 


19 


17 


3-y + p T + jets 


2 


1 









TABLE V: Number of signal events in the di-photon channel, 
after applying the basic cuts at an integrated luminosity of 
3 fb^ 1 and the center of mass energy of 7 TeV for all our 
benchmark points. 

The number of events in the tri-photon channel are 
relatively small since one of these photon comes from ra- 
diative decay of xi- The overall tri-photon event rates 
are small due to following two reasons: the smaller ra- 
diative decay branching fraction of the second lightest 
neutralino and together with the fact that the photon 
out of a X2 decay comes with relatively small px (see 
Fig. 5), because of small mass splitting between m^o 
and to^o . Hence in very small fraction of events they 
pass the requisite hardness cut. The effect is much more 
sever in case of four-photon channels as one can see from 
Table IV. Though we have quoted the event rates also 
in this case since one hardly has any contamination from 
the SM backgrounds, one still needs higher luminosity 
for better statistics. 



VII. SUMMARY AND CONCLUSION 

We have considered a supersymmetric scenario in 
which the gravitino (with a mass ~ 1 keV) is the LSP and 



the NLSP is the lightest neutralino. The second lightest 
neutralino is nearly degenerate in mass with the light- 
est neutralino. A possible origin of such a degeneracy at 
the low-scale lies in the form of non-universal high scale 
(~ 10 16 GeV) inputs of the soft SUSY breaking gaug- 
ino mass parameters. We have pointed out that such 
non-universal high-scale inputs can be realised in various 
representations of the SU(5), 50(10), and E(6) GUT 
group. 

We have examined the decays of the NLSP and the 
second lightest neutralino at the LHC. In such a scenario 
the second lightest neutralino has a substantial branching 
fraction of decaying into a photon and the lightest neu- 
tralino. The branching fraction depends on fj,, tan/3, and 
other scalar masses in the theory. The lightest neutralino 
is predominantly a bino and it too decays into a photon 
and a gravitino with a large branching ratio. Thus one 
naturally has spectacular multi-photon final states in a 
collider experiment, where light neutralinos are produced 
in abundance. The photons out of the NLSP decay are 
non-pointing and can be identified in the ATLAS inner- 
detector with an efficiency of 60%. Such non-pointing 
photons are free from any SM contamination. 

We have studied the di-photon, tri-photon, and four- 
photon final states in association with hard jets and miss- 
ing transverse energy in the context of LHC both at 
E cm = 7 TeV and 14 TeV and at an integrated luminos- 
ity of 3 fb _1 and 100 fb -1 , respectively. Though the di- 
photon and tri-photon signals look promising, one needs 
higher luminosity for the four-photon case. 

Detection of such multi-photon final states compris- 
ing non-pointing photons at the LHC would have serious 
implications for early- Universe cosmology and supersym- 
metry model building. On one hand one needs to have a 
suitable supcrsymmctry breaking mediation mechanism, 
which allows for light gravitino with a mass ~ 1 keV and 
nearly mass degenerate %° an d Xi- O n the other hand, 
this may give some hints towards a non-standard cos- 
mological scenario leading to a keV gravitino which is a 
warm dark matter candidate with right relic abundance. 
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